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Abstract 

We analyse the impact of growth conditions on asymmetric magnetic bubble expansion under 
in-plane field in ultrathin Pt / Co / Pt films. Specifically, using sputter deposition we vary the 
Ar pressure during the growth of the top Pt layer. This induces a large change in the interfacial 
structure as evidenced by a factor three change in the effective perpendicular magnetic anisotropy. 
Strikingly, a discrepancy between the current theory for domain-wall propagation based on a simple 
domain-wall energy density and our experimental results is found. This calls for further theoretical 
development of domain-wall creep under in-plane fields and varying structural asymmetry. 
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A. Introduction 


The recent manifestation of interfacial Dzyaloshinskii-Moriya interaction (DMI) [1-3] in 
nominally symmetric ultrathin Pt / Co / Pt films and Pt / [Co / NijAr multilayers [4-8] has 
raised questions concerning the origin of this interaction. In these multilayers DMI is usually 
explained in terms of structural inversion asymmetry, that arises as a result of the asymetric 
stacking of materials, hence, in symmetric Pt / Co / Pt the structural asymmetry should have 
a different origin. Most likely the effective DMI that is found in these structures arises from 
an asymmetry in the interfaces at either side of the ultrathin ferromagnetic layer [5, 9-12]. 
For Pt / Co / Pt samples even opposite signs of the DMI induced helds have been reported, 
where Je et al. [4] found an DMI induced held of fioHoMi ~ +26 mT, Hrabec et al. [5] found 
fioHoMi ~ -100 mT. Even more striking is that when replacing the top Pt layer by an Ir layer 
a reversal in H^mi is observed. This is counterintuitive as Pt / Co and Ir / Co interfaces 
are expected to have opposite DMI [8, 10], which should effectively increase the net Hdmi 
in a Pt / Co / Ir multilayer. The origin of these contradicting reports might lie in interfacial 
quality-dehning properties such as roughness, degree of intermixing, etc. Characterizing 
these structural properties quantitatively however, remains an outstanding challenge due to 
the ultrathin ferromagnetic layers used (typically < 1 nm) and poly-crystalline nature of the, 
typically sputter-deposited. Elms. Furthermore, in the area of magnetic held and current- 
induced domain-wall (DW) motion there have been widely dihering reports on the strength 
of the DMI [6-8, 13]. This might be due to intrinsic diherences between the deposited films 
grown in diherent laboratories. Moreover, any comparison is hampered since parameters 
such as growth rate, power, gas pressure and sample-substrate details are missing and, 
if given, are difficult to compare as the used specihc sputter apparatus and preparation 
method are also defining. Tackling these issues are key to understanding, comparing and 
interpreting the reported DMI results. The urgency of understanding the DMI in ultrathin 
hlms arises from the predicted huge impact on technological relevant devices where a chiral 
magnetization texture is preferred such as racetrack memories and logic devices based on 
DWs and skyrmions [7, 14-18]. 

Here, we investigate the ehect of growth conditions on the DMI by means of magnetic 
bubble expansion under in-plane helds in the archetype hlms of nominally symmetric Pt(4 
nm) / Co(0.6 nm) / Pt(4 nm) hlms as was successfully applied before [4, 5]. Specifically, us- 
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ing DC magnetron sputter deposition, we vary the Ar gas pressure during the growth of the 
top Pt layer ptop- Thereby, we change the growth kinetics of the top Pt layer [19, 20]. This 
leads to a different interfacial quality and/or degree of intermixing, and hence a variation 
in the degree of structural inversion asymmetry, giving rise to an effective interfacial DMI 
interaction [9], and magnetic properties. We hnd a large dependence of the perpendicular 
magnetic anisotropy (PMA) on ptopi varying up to a factor three between the lowest and 
highest ptop- Asymmetric bubble expansion under applied in-plane helds is observed indicat¬ 
ing a hnite DMI [4, 5]. We are, however, unable to describe our experimental data with the 
simple theory used before [4, 5]. This indicates a more complicated physical picture and the 
currently used models need to be expanded to incorporate the complex behavior observed. 
Our results shed light on the origin of the interfacial DMI in sputter deposited Pt / Co / 
Pt layers and provide a simple way to investigate the effect of changing interfacial quality. 

This letter is structured in the following order; we will start with introducing the used ex¬ 
perimental methods in section B. In section C we will present the basic magnetic properties 
of our samples. Here the large effect of the Ar pressure during the top Pt layer growth on the 
PMA is shown. In section D the results of expanding bubbles are presented, furthermore, 
we will conhrm that the DWs in all our samples follow the creep law. Furthermore, we will 
extract the creep law parameters and correlate this with the behavior of PMA as a function 
of Ptop- In section E we will concentrate on the results of bubble expansion under in-plane 
helds. We will elaborate on the differences and correspondence of the data with the current 
understanding. Finally, we will discuss the results and conclude in section F. 


B. Methods 

The samples are grown using parallel face-to-face target and substrate Ar DC magnetron 
sputter deposition from 2” targets in a system with a base pressure of ~ 3 x 10“® mbar. All 
samples are prepared on precut Si substrates (0.5x0.5 cm^) with a native oxide layer. The 
substrates were cleaned by acetone and isopropanol in an ultrasound bath followed by an 
in-situ 5W O 2 plasma exposure for 10 minutes prior to the deposition. The studied samples 
consist of SiO2//Ta(4)/PtL(4)/Co(0.6)/Ptu(4) (thickness between parenthesis is given in 
nm) where we have labeled the lower and upper Pt layers as PtL and Ptu. For Ta, Co and 
Pt DC magnetron powers of 20, 20, 60 W were used, respectively. The target-sample distance 
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during sputtering for Ta, Co and Pt was kept constant at 95, 95 and 195 mm, respectively. 
The Ar pressure during Ta, Co and PtL layer deposition was kept constant at 1.4, 1.0 and 
0.29 Pa, respectively. The pressure ptop during the Ptu layers was varied between 0.29 and 
2.8 Pa. The Pt growth rate was calibrated for every pressure to keep the top Pt layer 
thickness constant at 4 nm. The satnration magnetization was measured using a SQUID- 
VSM at room temperature. The perpendicular magnetic anisotropy was determined using 
angle dependent anomalous Hall effect measurements at room temperature. A standard 
polar MOKE setup was used to obtain the hysteresis loops at a constant field sweep rate 
of 10 mT/s. The bubble expansion measnrements were performed using wide held Kerr 
microscopy in polar mode with a custom perpendicular pulse coil with a calibrated rise time 
of 70 jUS to provide the perpendicular to the surface z-£elds and a standard in-plane magnet 
for the in-plane x-held. 


C. Basic magnetic properties as a function of ptop 

In Fig. 1(a) the measured saturation magnetization Ms as a fnnction ptop is plotted for 
Si02 // Ta(4 nm) / Pt(4 nm) / Co(0.6 nm) / Pt(4 nm). The Ms is found to be close to the 
bnlk value of Co (1440 kA/m) for all used ptop- In Fig. 1(b) we plot the effective anisotropy 
held HK,eff as a fnnction of ptop which shows a clear increase with ptop up to a factor 3.1 
between ptop = 0.29 Pa (760±80 mT) and ptop = 2.80 Pa (2380±240 mT). 

In Fig. 1(c) the hysteresis loops of the samples are presented showing sqnare loops and 
increasing coercivity for higher ptop, in agreement with observed behavior of Hx,eff- 
These results show a clear trend of the perpendicular magnetic anisotropy (PMA) as cap- 
tnred in i7x,e// with increasing ptop illustrating the sensitive nature of the PMA to the 
growth conditions. Hence, the structural quality of the Co layer and interfaces with the Pt 
are paramonnt and are determined by the kinetics dnring the growth of the top Pt layer 
[19, 20]. Similar behavior has been seen before where the PMA was fully dominated by the 
bottom Pt / Co interface [5]. In another report the lack of PMA from the top Co / Pt was 
attributed to a strongly intermixed top Co / Pt interface [21]. Hence, extending on this 
scenario a basic nnderstanding can be found by examining the incoming kinetic energy of 
the Pt atoms dnring growth, which ronghly scales as the inverse of the pressure [19, 20]. 
To test this we plot in Fig. 1(d) HK,eff as a fnnction of Indeed for high pressnre (> 
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0.8 Pa, < 1.19 Pa“^) a linear scaling with is found as indicated by the linear £t (red 
dashed line). For lower ptop this trend breaks down. We attribute this break down as the 
region where there is a strong intermixing between the Pt and Co atoms possibly creating 
an alloyed PtCo layer, which are known to show high PMA. In this case no well-defined Pt 
/ Co and Co / Pt interfaces are formed, i.e. an interfacial structural inversion asymmetry 
would be poorly defined. 

For higher pressures there is less intermixing and the interfaces get increasingly better 
defined. In this regime an effective DMI due to interfacial structural inversion asymmetry 
could be expected as long as the bottom and top interface are sufficiently different. We have 
labeled these regions the ‘alloying’ and ‘layered growth’ regime, respectively, in Fig. 1(d). 
Other origins for the increasing PMA with ptop could be a changing density of the top Pt 
him with Ptop, different crystal orientations, growth modes and/or grain sizes. 

We have performed cross-sectional and plain view high resolution transmission electron 
microscopy (HRTEM) studies on samples prepared on SiN^. membranes and 50 nm thick 
lamella prepared using a focused ion beam but were unable to discern differences between 
the samples grown with low and high ptop, both in crystal size as in orientation, viz. a strong 
<111> texture for all samples was observed from electron diffraction patterns. Moreover, 
using scanning TEM - EDS (energy dispersive X-ray Spectroscopy) mapping we were unable 
to discern any differences in the intermixing between Pt and Co. This is due to the limit 
in resolution in the element maps due to the combination of the limited thickness of the 
ultrathin Co him (0.6 nm) and the surface roughness of the underlying Pt layer (due to 
its poly-crystalline nature, yielding a projected layer thickness of the Co layer of ~1.5 nm. 
Future atom-probe tomography measurements might shed light on the structural properties 
and provide further proof for the given hypothesis. We conclude that increasing the Ar 
pressure during the growth of the top Pt layer in Pt / Co / Pt layers increases the effective 
PMA by a factor 3.1. We attribute this change to an increasingly better defined layered 
growth of the Pt / Co / Pt stack. 


D. DW creep - bubble expansion 

In Fig. 2(a) a typical example of a symmetrically expanding bubble is shown under 
application of held pulses using the following measurement sequence (identical expansion 
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is observed for all samples). (1) The sample was saturated with a z-£eld Hz^s pulse with 
a length of ts- (2) A background image is taken and subtracted. (3) A circular bubble is 
then nucleated using a nucleation pulse {Hz^N,tN)- (4) The in-plane field was set. (5) 
An initial reference image is taken. (6a) A propagating pulse {Hz^p,tp) is applied. (6b) 
An image is recorded. Step 6 is then repeated X times. (7) is switched off. The 
obtained images are then processed to extract the bubble nucleation origin and DW velocity 
along the bubble circumference. Special care was taken to find a nucleation site where 
a consistent and circular bubble was nucleated with both polarities during the nncleation 
pulse. Furthermore, the area around the nucleation site, i.e. the area where the expanding 
bnbble travels throngh, was checked for strong pinning sites which conld possibly lead to 
non-continnous extracted DW velocities. After monnting every sample the in-plane held 
coil angle relative to the sample plane was carefnlly tnned to minimize the leak into the 
perpendicular held direction. 

By fitting and averaging the expansion of the bnbble versus the total applied Hz pulse 
time the DW velocity v versus iiqHz is extracted, as shown in Fig. 2(b) (black symbols). 
The data can be fitted to the creep law [22]: 

n = noexp(-x(/ioi72)"^/'^). (1) 

The red dashed lines in Fig. 2(b) are fits to the creep law and an excellent agreement is 
fonnd. The bine symbols correspond to the top and right axis where we have plotted ln(n) 
versus (/iohfz)”^^^ showing a linear behavior, as expected. Identical behavior is observed in 
all samples as is shown in Fig. 2(c). In Fig. 2(d) we plot the characteristic fit parameters 
uo (left inset) and y as a fnnction of ptop obtained from fits to the data shown in Fig. 2(b). 
Fnrthermore, Uq is the characteristic speed parameter and scales with In fhe creep 

description vq is defined as the prodnct of the correlation length ^ of the disorder potential 
(average distance between pinning sites) and the depinning freqnency /o [23]. From this 
data, however, we cannot distinguish which of these two parameters changes as a fnnction 
of Ptop- Speculative, we expect the change in ^ with ptop to be small as we observe identical 
structural properties for all ptop from the HRTEM study and the strong increase in vq could 
be dominated by /o, this however, remains an open issne. To fnlly nnderstand this relation 
further theoretical investigation is needed. 

The scaling constant y, scales with Hx^eff which is plotted in red in Fig. 2(d) (right 
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axis). X is originally defined as UcH'^^^^/{ksT) in the DW creep theory with /i = 1/4 [22, 24], 
where Uc is an energy scaling constant, is the critical magnetic depinning held at zero 
temperature, and ksT denotes the thermal energy. Using the description of the creep law 
X ~ -^Xe// obtained [25]. In the right inset of Fig. 2(d) we have plotted x versus 

HK,eff and the dashed red line shows x ~ -^xe//- good agreement leads us to conclude 
that no large changes in the shape of the DW pinning landscape is expected, as it only scales 
with Hx,eff- Summarizing, we conclude that the DW motion is well characterized by the 
DW creep law for all ptop and Hz, i.e. velocities, studied. 


E. Bubble expansion under in-plane field 

In Pt / Co / Pt thin dims with PMA, Bloch type DWs are preferred due to the magneto¬ 
static DW anisotropy held Hk 4 w = 4:Kdw/i'n'PoMs) with Kdw = NxfJ,oMg/2 the magneto¬ 
static DW anisotropy, where W = tco ^og{2)/{nA) is the demagnetization coefficient of 
the DW [26] with tco the Co layer thickness and Ms the saturation magnetization. For 
magnetic bubbles this leads to the situation as shown in the top left panel of Fig. 2(e) 
where a magnetic bubble with its magnetization along +z (grey area) is shown. Since 
there is no preferred chirality for Bloch DWs the DW magnetization can rotate either clock 
or anti-clockwise going through the DW from the inside to the outside of the bubble as 
indicated by the double arrows. In a simple description of DMI at DWs [4, 5, 10, 13], the 
DMI manifests itself as a built-in magnetic held Hr ,mi = D/[p^MsA) pointing in-plane 
perpendicular to the DW, where D is the strength of the DMI interaction, and A the DW 
width. Hence, for Hbmi > Hk^w Neel type DWs are preferred as shown in the top right 
panel of Fig. 2(e) and depending on the sign of D a clock (D > 0), or anti-clockwise {D < 0) 
DW chirality is introduced as indicated by the white and purple arrows, respectively. For 
Hdmi < Hk 4 w the DW assumes a mixed Bloch-Neel character. By now applying a strong 
in-plane held along the x-axis (Ffx) the DWs magnetization reorients itself along due to 
the Zeeman energy. This is shown for a bubble with D > 0 in the bottom left panel of Fig. 
2(e), where we assume H^mi > Hk 4 w The DW segments that are parallel to (top and 
bottom of bubble, orange arrows) become Bloch type DWs with an increase in DW energy 
density due the Zeeman energy. The DW segments that are perpendicular (left and right of 
bubble) remain of the Neel type. The Neel DW which has its direction reversed (and thus its 
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chirality) however, undergoes an increase in DW energy density. The same happens when 
D > 0 (bottom right panel of Fig. 2(e)), albeit here the left and right DW behavior are 
reversed. Hence, an in-plane magnetic held breaks the symmetry of the DW energy prohle 
along the in-plane held direction. 

This is experimentally shown in Fig. 2(f) where the bubble expansion for the ptop = 1.40 
Pa sample is shown with a held applied. A strongly asymmetric expansion is observed 
where the DW moving in the direction of the applied held moves mnch faster as the 
DW moving against the in-plane held. The inset shows the expansion in the same sample 
bnt with inverted Hz showing identical asymmetric expansion albeit mirrored in the y- 
axis. This observation was attributed to a built-in DMI held manifesting itself at the DWs 
as explained before [4, 5]. By applying an in-plane held dnring the expansion the radial 
symmetry is broken as the DMI held would prefer Neel type DWs with a certain chirality 
which is broken by the in-plane held. 

In Fig. 3(a)-(f) we plot the obtained DW velocities n as a fnnction of along the 
extremes of the bnbble along the x and y-axis, i.e. parallel and perpendicular to the applied 
\HJ^ held for all ptop- Note the diherent v scales used; the variation in v increases with 
higher ptop. We have averaged the left moving DW segment velocity with positive > 0 
(i.e. np-down DWs when moving from inside the bnbble to ontside), which we label L-|-, 
with the right moving DW velocity for < 0 (R-) (i.e. down-np DWs when moving from 
inside the bnbble to outside), as these have to be identical by symmetry (see Fig.2(f)). 

The same has been done for the np (U) and down (D) moving DW segments. This allows 
us to compensate for small misalignments of the, relative to Hz, large Hx held. Despite this 
compensation small residual asymmetries remain which we consider negligible (compare the 
U-l-D- (green) and U-D-l- (blue) prohles which need to be similar by symmetry). Extracting 
the DW velocity prohles in this way allows us to disentangle the in-plane held ehect on the 
DW velocity parallel (L-l-R- and L-R-I-) and perpendicular (U-l-D- and U-D-I-) to the applied 
Hx held. The asymmetric bnbble expansion is rehected in the diherence between the L-l-R- 
and L-R-l- velocity prohles at a certain Hx- The U-l-D- and U-D-l- prohles show a symmetric 
prohle aronnd Hx = 0 mT. This observation corresponds with the aforementioned picture 
of a DMI held breaking the in-plane symmetry of the bnbble expansion. 

For Ptop = 0.29 Pa an increasingly ill-dehned bnbble expansion and extra domain nu- 
cleation is observed for high \Hx\ (not shown). This makes it impossible to extract well 



defined DW velocities for \Hx\ > 200 mT as there is no bubble expansion but a dendrite like 
magnetization reversal. On close inspection of all the prohles three distinct DW velocity 
prohle characteristics can be identified: (1) symmetric bubble expansion perpendicular to 
the in-plane held direction (U-l-D- and U-D-I-). (2) Asymmetric bubble expansion parallel 
to the in-plane held direction (L-l-R- and L-R-I-). (3) For = 0.29 Pa the DW velocity 
shows a velocity drop symmetrically around Hx = 0 (most clearly seen for the U-l-D- and 
U-D-l- prohles) which decreases for the samples grown with higher ptop- 

All the experimental L-l-R- (black symbols) prohles show a minimum at ~ 60 mT. In 
the interpretation of an ehective DMI held this would indicate a constant positive built-in 
DMI held Hdmi > 0, i.e. a positive DMI energy D > 0 would be found for all samples 
preferring a right handed DW chirality. Furthermore, there would be no variation in the 
DMI as a function of ptop which is puzzeling as the large variation in PMA indicates a strong 
change in the interfacial quality and speculative also on the structural inversion asymmetry. 
The found D > 0 is in-line with the report of Je et al. but opposite to the report of Hrabec 
et al. for Pt / Co / Pt. Moreover, the overall v prohle shape observed here is very diherent 
than described by the theory used in these reports. This leads us to conclude that the 
interpretation of the DMI leading to a simple in-plane held might be too simple for our 
samples. 

The out-of-plane driving held might have an ehect on the prohle shapes as this determines 
the strength of the DW motion. To examine the ehect of the out-of-plane driving held we 
have plotted the Inn L-l-R- prohles in Fig. 3(g)-(l) for diherent \Hz\ drive helds. For the 
lowest Ptop we observe an overall changing shape of the velocity prohle with increasing \Hz\, 
for the higher ptop the overall shape remains constant. This behavior further hints towards 
diherent mechanisms at work between the low and high ptop samples which might be due 
to diherent DMI strengths. From this data we can extract the no and y parameters by 
htting the velocity dependence on \Hz\ to the creep law (Eq. 1) for every Hx- This is shown 
in Fig. 3(m)-(r) where we plot In(no) and y as a function of Hx- The cyan line in the 
In(no) data is a ht to In(no) = a -|- b\Hx\, which shows that a symmetric behavior for vq is 
observed relative to Hx = 0 for all samples. Furthermore, the variation in In(no) decreases 
strongly with increasing ptop and remains more or less constant for the highest ptop, which 
indicates a changing underlying dynamic response captured in uq- The y data, however, 
shows an increasing asymmetry relative to Hx = 0 with higher ptop. The origin of this effect 
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is not clear at the moment but the overal shape of the x{Hx) dependence is very different as 
expected from the theory described by Je et ai, where a simple inverted parabolic behavior 
is expected with a maximum in y at —Hjjmi- 


F. Discussion & Conclusion 

Due to the discrepancy between the previously used theory and our experimental data we 
are unable to fully interpret the experimentally obtained velocity prohles. Naively interpret¬ 
ing the held at which a minimum in the DW velocity was observed, we hud ~ 

60±10 mT for all ptop- If we interpret the interfacial Hdmi as an effective held, i.e. origi¬ 
nating from the diherence of D between the bottom Pt / Co and top Co / Pt interfaces, a 
correlation between the PMA (through ptop) and Homi should have been evident. This is 
motivated by the known dependence of the PMA on the interface quality. As this is not the 
case, i.e. a constant H^mi is found for all ptop whilst the PMA varies over a factor three, 
indicates that this interpretation might be hawed. Hence, this calls for further theoretical 
development of the underlying mechanisms and the simple model introduced by Je et al. 
should be used with great care. Specihcally, we suggest that full micro-magnetics of the 
DWs prohles and the strong pinning of the DWs in the creep regime should be considered. 
The successful application of the theory in the reports by Je et a/.[4] and Hrabec et al.[5] 
might be due to diherent preparation procedures and (lower) held regimes studied. The 
Co layer in our samples are rather thick (0.6 nm) compared to the sample used by Je et 
al.[4] (0.3 nm) which exhibit ultra low DW pinning, but are similar to the samples grown 
by Hrabec et al.[5]. 

Very recently, asymmetric bubble expansion under in-plane helds was attributed to a 
completely diherent mechanism, viz. chiral-dependent damping [27] again in similar sys¬ 
tems of Pt / Co / Pt. This would modulate the attempt frequency Vq in the creep law (Eq. 
1) as a function of in-plane held and was explained as a dissipative (held-like) spin-orbit 
torque on the DW dynamics. As we have seen from the data in Fig. 3(m)-(r) we observe 
a symmetric behavior of vq relative to = 0 and can exclude this in our samples. In 
fact, the experimental hnding that y, and not Inuo, is asymmetric with respect to in-plane 
held suggests that the main asymmetry is in the DW energy rather than its dissipation. 
This conclusion holds provided the collective-pinning theory of the creep regime applies 
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but is otherwise model-independent. One possibility for the influence of in-plane fields on 
the DW energy not considered before is the tilting of the domains by the in-plane held. 
This might be signihcant as the in-plane helds we apply are not necessarily small compared 
to the PMA. The tilting of the domains by the in-plane held leads to an ehective reduc¬ 
tion of the driving force on the DW. In the simplest model for this reduction one replaces 

in the creep formula. Using this model with HK,eff as a ht 
parameter, we obtained htted values for Hx^eff fhat have the same order of magnitude as 
the ones that we measured directly, but were nonetheless unable to accurately describe all 
our measurements. This suggests that a more accurate calculation of the DW energy, in¬ 
cluding the tilting of domains by the in-plane held, is needed, which is the subject of future 
theoretical work. 

In summary, we studied the DW velocity prohles of expanding bubbles in diherently 
prepared Pt / Co / Pt samples under the application of in-plane helds. Specihcally, we 
varied the Ar pressure during the deposition of the top Pt layer leading to a factor three 
increase in PMA between the lowest and highest pressure used. This indicates a large change 
in the structural quality between the samples. The velocity prohles of the expanding bubbels 
were successfully extracted. However, the results could not be described by a model that was 
used successfully before. We believe that the found results shed light on the understanding 
of the ehective interfacial DMI in ultrathin hlms and will facilitate the quest to boost the 
DMI to stabilize chiral Neel DWs and create skyrmionic spin textures in thin him systems 
at room temperature. 
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FIG. 1. (a) Saturation magnetization Ms as a function of ptop- (b) Effective anisotropy field 
poHK,eff as function of ptop- (c) Polar MOKE loops for different ptop- (d) Effective anisotropy 
field pQHx^eff as a function of the red dash line is a linear fit to the data for p^^ <1.19 Pa“^. 
The blue line indicates a transitions between the ‘layered growth’ and ‘alloying’ regime. 
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p,„p = 0.29 Pa 

H=0rr\TH,p = + 5.39 mT, t. = 9.93 ms 







FIG. 2. (a) Typical example of a symmetrically expanding bubble with = 0 mT driven by an 

out-of-plane magnetic field of 5.39 mT with a duration of tp = 9.93 ms for the ptop = 0.29 Pa 

sample, (b) Typical velocity v versus PqHz, black symbols and left-bottom axis and ln(u) versus 

blue symbols and right-top axis for ptop = 0.29 Pa. The red dashed lines are a fit to the 

creep law. (c) ln(u) versus {poHz)~^^^ for all samples, (d) Creep law scaling constant y versus ptop 

(black symbols and right/bottom axis). Effective anisotropy field as function of ptop (red symbols 

and right/bottom axis). The left inset shows the creep law prefactor vq as a function of ptop- The 

right inset plot x versus Kf,f j indicating the x ~ -^e// Proportionality, (e) Top left panel shows a 

Bloch bubble where we have indicated the magnetization angle inside the DW defining the bubble. 

The top right panel shows a Neel type bubble where the chirality of the DW depends on the sign 

of the DMI interaction. Bottom two panels show the orientation when the DWs magnetization is 

saturated along the applied in-plane field direction and the thickness of the arrows indicate the 

increase in energy of the DWs for different D, the white dotted circles in the bottom right panel 

accentuate the change when the in-plane field is reversed, (f) Bubble expansion (same scaling as 

in (a)) in the ptop = 1-40 Pa grown sample while applying an in-plane field, the inset shows the 
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expansion with inverted Hz showing identical, but mirrored in the y-axis, expansion. 
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FIG. 3. (a - f) DW velocity profiles for L+R-, L-R+, U+D- and U-D+ on a log scale as a function 
of Hx for all ptop- Note the different ranges of the v scale used, i.e. the velocity variation increases 
with increasing pressure ptop- (g)-(l) L+R- velocity profiles as a function of Hx for all ptop and 
different \Hz\ as indicated in the panels, (m)-(r) In(uo) and x as a function of Hx for all ptop as 
extracted from the data shown in (g)-(l). The cyan line in the Inuo graph indicates the symmetric 
behavior relative to Hx = 0. 
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